On the critical behavior of two-dimensional liquid crystals 
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The Lebwohl-Lasher (LL) model is the traditional model used to describe the nematic-isotropic 
transition of real liquid crystals. In this paper, we develop a numerical study of the temperature 
behaviour and of finite-size scaling of the two-dimensional (2D) LL-model. We discuss two possible 
scenarios. In the first one, the 2D LL-model presents a phase transition similar to the topological 
transition appearing in the 2D XY-model. In the second one, the 2D LL-model does not exhibit 
any critical transition, but its low temperature behaviour is rather characterized by a crossover 
from a disordered phase to an ordered phase at zero temperature. We realize and discuss various 
comparisons with the 2D XY-model and the 2D Heisenberg model. Adding to previous studies 
of finite-size scaling behaviour of the order parameter and conformal mapping of order parameter 
profile, we analyze the critical scaling of the probability distribution function, hyperscaling relations 
and stiffness order parameter and conclude that the second scenario (no critical transition) is the 
most plausible. 

PACS numbers: 64.70.M-, 64.60.Bd, 64.70.rQf, 05.70.Jk, 05.50.+q, 68.35. Rh 



I. INTRODUCTION 

At high temperatures, the natural state of a liquid 
crystal is the isotropic phase, because all the molecule 
orientations are equally probable. At lo-w temperatures, 
an order along a preferred direction may appear, it is 
the nematic phase. When such a material is cooled at 
even lo-wer temperatures, other types of ordered phases 
are likely to appear, the description of -which requires re- 
alistic potentials (see e.g. reference jlj]). The transition 
bet-ween the isotropic and the nematic phases generally 
occurs at a critical temperature. Various experiments on 
three-dimensional liquid crystals, exhibited a -weak first 
order transition 2]. Ho-wever, the corresponding two- 
dimensional (2D) problem is far from begin that clear. 
Experimentally, only 2D liquid crystals composed of long 
rigid interacting rods have been studied. 

First, the 2D problem must be defined through the 
symmetries: there are generally a continuous global sym- 
metry 0{n), and a local Z2 symmetry. For example, the 
t-wo most-studied models deal -with 0{2) (the XY-model), 
or 0(3) (the Heisenberg model) symmetries. The Hamil- 
tonian reads in both cases as: 



n = 



Sj, 



(1) 



-with S,; a n-components spin, of unit length, at the lat- 
tice site i. The sum runs over all pairs of neighbouring 
spins of the 2D square lattice and J > is the spin-spin 
coupling constant. 

According to the Mermin-Wagner-Hohenberg theorem 
d, 0, H, H, Q J the order parameter fluctuations prevent 



from any spontaneous breakdo-wn of continuous symme- 
tries for 2D systems -with short-range interactions. Ho-w- 
ever, this theorem cannot conclude in the case of transi- 
tions to a phase with no long range order. An example of 
such a 2D critical behaviour, which is not a continuous- 
symmetry breaking, is the 2D XY-model {0{2), Abelian 
symmetry). This model presents a topological transi- 
tion characterized by a critical low temperature phase 
as a plasma of vortices pairs, at temperatures under 
the Berezinskii-Kosterlitz-Thouless transition tempera- 
ture, Tbkt [1, 0- At much lower temperatures, the 
system state is dominated by spin waves. The high tem- 
perature phase occurs when the pairs dissociate. 

On the other hand, the 2D Heisenberg model (0(3), 
non-Abelian symmetry) has no transition at any finite 
temperature (asymptotic freedom) [Tol , 

El, The 

longitudinal modes of the 0(3) model are frozen for the 
low temperatures, and only the transverse modes are ac- 
tivated, leading essentially to two Gaussian modes. The 
deep difference between the two models at finite tem- 
peratures, is the stability of the vortices. Indeed, the 
topological transition of the 2D XY-model is the result 
of the stability of the topological defects, while the "third 
spin dimension" in the 0(3) model, makes the vortices 
unstable at any finite temperature. 

However, the above picture has been questioned re- 
cently by the numerical evidences for a possible transition 
in the 2D Heisenberg model [M [H, [ll, [13] ■ Moreover it 
has been shown that a quasi long range order (QLRO) 
phase [isl . IToj might appear at very low temperatures in 
finite systems. Indeed in these numerical works, the ther- 
modynamic limit remains questionable and the analytic 
approach of [l^ relies explicitly on the finite size of the 
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lattice. 

A model for a regular 2D liquid crystal system was pro- 
posed by Lebwohl and Lasher (LL) [2^. Based on a lat- 
tice version of the mean field theory of Maier and Saupe 
[2l| . the molecules are represented by n— component unit 
vectors and the "spin-spin" interactions are given by the 
second Legendre polynomial, P2 (keeping the local Z2 
symmetry) . The popularity of the LL- model comes from 
its ability to reproduce the weak first order phase transi- 
tion observed experiinentally between the isotropic and 
the nematic phases [2]] in the three-dimensional space. 
Then, the question of possible phase transitions in this 
model is attractive and has been addressed in numerous 
studies. ;i8, 22, 2|, [H, [H, [H 

In the past twenty years, several numerical evidences 
pointed out a possible topological transition for the 2D 
nematic-isotropic phase transition in liquid crystals. The 
situation is clear in the case of planar rotator models 
(see e.g. references [13, HI]), but still controversial for 
three-component models. In reference [l^l, we have stud- 
ied the nematic-isotropic transition of 2D liquid crys- 
tals using a 0(2) vector model characterized by nonlin- 
ear nearest-neighbour spin interactions governed by the 
fourth Legendre polynomial P4. The system has been 
studied through standard finite-size scaling and confor- 
mal rescaling of the density profiles. We also estimated 
the Binder cumulant [2^ as a function of the tempera- 
ture for different values of system size L. This cumulant 
has been proved to be universal at a critical tempera- 
ture. Evidences for a topological transition at a finite 
temperature have been underlined. 

In the case of the three-component model, in a remark- 
able work Kunz and Zumbach [2^ concluded in favour 
of such a critical topological transition. Although they 
performed a careful and sizable study, they were unable 
to decide conclusively between essential singularities or 
standard power laws for the correlation length and the 
susceptibility, when approaching the critical temperature 
from the high temperature phase. However, they devel- 
oped a qualitative picture (pairing of topological defects 
which carry most part of the system energy, sharp in- 
crease of the defects density and seeming discontinuity of 
the rotational-rigidity modulus, finite cusp of the specific 
heat, proliferation of unbounded defects at high tempera- 
ture) and their conclusion inclined clearly to the essential 
singularity, for a temperature Tbkt- More recently, we 
studied the 2D LL- model [H, [2^, HI] , and found signals 
in favour of a QLRO phase at T < TbkTi with a magne- 
tization that decays with the system size as a power law 
with the critical exponent ri{T)/2 26]. We obtained good 
agreement with the ri{T) exponent obtained using the 
powerful technique of the conformal transformation [25^ . 
In [1^, the value of 77 (T) was estimated using directly 
the susceptibility, and the agreement appeared to be rea- 
sonable. Possible discrepancy could come from the finite 
size of the systems. This QLRO phase could be similar to 
the low temperature phase of the 2D XY-model. Mondal 
and Roy [SO], using standard finite size scaling method 



and Monte Carlo simulations for the 2D LL-model, gave 
arguments for a continuous phase transition. The criti- 
cal temperature was estimated to be Tc = 0.548 ± 0.02, 
a httle higher than 0.513 calculated in [H, [H, [H [11] . 
A still more recent study reported evidences in favour of 
a topological transition '31^. Then, the transition might 
be driven by stable topologically point defects known as 
^-disclination points. 

Despite 20 years of numerical results claiming for a 
topological transition for the 2D LL-model, a recent work 
called seriously this scenario into question. Indeed, pre- 
cise estimates of the Binder cumulant for various sys- 
tem sizes, have been unable to show any crossing point 
at a definite temperature [12] • This result contradicts 
completely the transition scenario at a critical point. In 
the same work 132], the possibility of a QLRO phase at 
T — 0.4 (a value well below the "transition" tempera- 
ture estimated previously for this model) was questioned. 
Two strong evidences were given to conclude that the 
2D Lebwohl-Lasher model does not show any quasi-long- 
range ordered phase. The evidences were the violation 
of the hyperscaling relation between the apparent expo- 
nents of the magnetization and the susceptibility, and the 
failure of the first-scaling collapse [l^ of the probability 
distribution function of the order parameter. 

The main goal of the present work is to present a com- 
plementary investigation of the 2D LL-model. We will 
compare with the same tools this model with the 2D 
XY- and Heisenberg models, for which the critical be- 
haviours are not questionable. Since the 2D LL-model 
seems to share some common properties with both the 
XY- and Heisenberg models, it appears helpful to pro- 
vide the comparisons. In section IIIII we will analyze the 
system with and without an applied magnetic field using 
the finite size scaling method (FSS). The results which 
follow from conformal transformation method (CT) will 
also be reminded [25j . We shall recover the apparent ev- 
idences for a sort of topological transition similar to the 
Berezinskii-Kosterlitz-Thouless transition. However, we 
will see in the section |Vl with new results on the shape 
of the tail of the probability distribution function for the 
magnetization, check of the hyperscaling relation, and 
study of the stiffness, how improbable is the reality of 
a critical topological transition, confirming the results of 



II. DEFINITIONS OF THE MODELS 

We will consider hereafter the 2D XY-model [1], the 
2D Heisenberg model [l^], and the 2D Lebwohl-Lasher 
model 20]. The simulations will be performed using 
Monte Carlo (MC) Wolff cluster algorithm ^34j appro- 
priately adapted to the model under consideration [22| . 
in order to avoid most of the critical slowing down be- 
haviour close to transition, if any. All simulations are 
realized on a square lattice of size L x L with periodic 
boundary conditions. The system sizes and number of 
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MC iterations will be specified when needed in the figure 
captions. The constants J and fc^ are fixed to unity. 

• For the 2D XY-model, the N = L"^ classical 
spins are confined in the lattice plane x, y. The 
spin Si is parameterized by: {Sf — cos 6i,Sf — 
sin^i), with 0i the angle between the direction 
of the spin Si and the x-axis. According to ((1]), 
the system state is governed by the Hamiltonian 
H = —J^^ij)C0s{9i — 9j), with the ferromag- 
netic coupling constant J > 0, and the sum run- 
ning over all the nearest-neighbour pairs of spins. 
The critical features are eventually characterized 
by the power-law behaviour of the magnetization 
per site: m = jf\/{J2i which is a non- 
negative real number. Thermal fluctuations of the 
order parameter give access to the susceptibility 
X = i7fcBT((m2)-(m)2). 

• For the Heisenberg model, the N = L"^ classical 
spins can point in any direction of the three di- 
mensional space x,y,z, (Sf = sinOi cos (pi, Sf = 
sin0ismipi, Si = cos9i), with the standard defi- 
nitions for the local angles. Similar to the above 
XY-model, the Hamiltonian has the form H = 
—J^(^ij^ Si • Sj, but with three-component spins 
and the local order parameter is defined by the ap- 
propriate magnetization per site. 

• For the Lebwohl-Lasher model, the N = classi- 
cal spins can also point in any of the three dimen- 
sional space directions x, y, z. The Hamiltonian of 
the system, 

W = -J^P2(S.-S,), (2) 

is an obvious generalization of the Heisenberg inter- 
action which guarantees the local Z2 symmetry of 
liquid crystals. P2 is the second Legendre polyno- 
mial, that is: ^2(2;) = (Bx^ — l)/2. One defines the 
local order parameter (also called nematization) as: 
m = ;^ ^2(cosai), where is the angle be- 
tween the direction of the spin S^ and the direction 
of symmetry breaking. 

III. QUALITATIVE ANALYSIS OF THE ORDER 
PARAMETER AND THE SUSCEPTIBILITY 

In the following we study the behaviour of the order 
parameter (magnetization or nematization) m and the 
magnetic susceptibility (or response function) Xi as a 
function of the temperature for various system sizes L. 

A. The 2D XY-model 

In figure m , the magnetization m and the susceptibil- 
ity X (top and bottom) for the 2D XY-model are plot- 



ted versus the temperature. Both plots have the same 
temperature range, wide enough to cover the domain of 
temperatures where important thermodynamic changes 
are expected to occur (essentially: the BKT transition, 
which should appear at the temperature Tbkt — 0.893). 
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FIG. 1: Temperature dependence of thermodynamic quanti- 
ties for the 2D XY-model. Top: magnetization (m). Bot- 
tom: susceptibility (x). (Monte Carlo simulations with 10^ 
MCS/spin, after removing the first lO" MCS to achieve proper 
thermalization) . 

We can observe in figure ([T]) a progressive decrease of 
the magnetization with the system size, at constant tem- 
peratures. Moreover, as become evident from figure ([2|. 
we see that these decays follow power-laws with the sys- 
tem size in the whole range T < Tbkt, according to the 
FSS picture in the critical phase, m L~2''(^), with an 
exponent r]{T) which depends on the value of the tem- 
perature (to be discussed in the next section) . When the 
value of the temperature is above Tbkt, the decay of the 
magnetization is faster than a power law. 

On the other hand, the magnetic susceptibility at fixed 
temperature (figure [T]) exhibits regular increase with the 
system size for all temperatures at and below Tbkt- It 
is in contrast with an ordinary second order phase transi- 
tion in which the susceptibility diverges only at the criti- 
cal temperature. In particular, the susceptibility exhibits 
power law behaviours with L for any T < Tbkt (figure 
[3]). It is interesting to notice that a change in the ex- 
ponent of the power-law of the susceptibility with the 
system size, appears at Tbkt- 

Also, we see in figure (IT|), that the susceptibility 
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FIG. 2; Behaviour of the magnetization versus the system 
size at different temperatures for the 2D XY-model. Power- 
law behaviours appear in the whole range T < Tbkt- (Monte 
Carlo simulations with 10® MCS/spin after removing the first 
lO" MCS to achieve proper thermalization) . 




FIG. 3: Behaviour of the magnetic susceptibility versus the 
system size at various temperatures for the 2D XY-model. 
(Monte Carlo simulations with 10^ MCS/spin after removing 
the first 10® MCS to achieve proper thermalization). 



reaches its maximum value systematically above Tbkt- 
The temperature at which this maximum occurs tends to 
Tbkt when L — > oo. 

The power law behaviours of the magnetization and of 
the magnetic susceptibility at low temperatures, can be 
interpreted as evidences of a QLRO phase for the 2D XY- 
model in the thermodynamic limit, in the temperature 
range T < Tbkt- 



B. The 2D Heisenberg model 



The 2D Heisenbergmodel is expected not to exhibit 
any phase transition jlOL [ill [l^ hj. The magnetic sus- 
ceptibility strongly diverges when the temperature van- 



ishes However, some authors argued in favour of a 
transition Il6l Il7l | and possibly an effective quasi- 

long-range-ordered phase [H, [l^ at very low tempera- 
tures, at least in large but finite systems. 
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FIG. 4: Temperature dependence of thermodynamic quanti- 
ties for the 2D Heisenberg model. Top: magnetization (m). 
Bottom: susceptibility (x)- (Monte Carlo simulations with 
10® MCS/spin after removing the first 10® MCS to achieve 
proper thermalization). 

In figure ([4]) , one can see the behaviours of the magne- 
tization and of the magnetic susceptibility with the sys- 
tem size and temperatures for the 2D Heisenberg model. 
We can notice the decrease of the values of the inflection 
point of TO versus T, when the system size increases. On 
the same way, the temperatures where the susceptibility 
reaches the maximum, are displaced to lower and lower 
values when the system sizes are larger. As a matter 
of fact, if there is a positive critical temperature analo- 
gous to the Tbkt, its value should be quite small, at the 
thermodynamic limit. 



C. The 2D LL-model 

We analyze now the magnetization and the magnetic 
susceptibility of the 2D Lebwohl-Lasher model, using the 
same procedure as for the 2D XY-model. In figure ([5]) 
the magnetization, to, and the magnetic susceptibility, x, 
are shown versus the temperature for various lattice sizes 
L. There is a complete qualitative agreement with the 
behaviour of the same quantities for the 2D XY-model 
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(figure [T]). Indeed the "critical" temperature is shifted to 
T* = 0.513 [H, d^, but all the power-law behaviours are 
equally recovered (figures [6] and [7]), and even the change, 
at r*, of the apparent value of the exponent of the sus- 
ceptibility versus the system size, or the asymptotic con- 
vergence to T*, of the temperatures at which the sus- 
ceptibility reaches its maximum are identically noticed. 
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FIG. 5; Temperature behaviour of thermodynamics quanti- 
ties for the 2D Lebwohl-Lasher model. Top: magnetization 
(m). Bottom: magnetic susceptibility (x). (Monte Carlo sim- 
ulation with 10® equilibrium steps (measured as the number 
of flipped Wolff clusters) after removing the first 10® MCS to 
achieve proper thermalization) . 
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FIG. 6: Behaviour of the magnetization with the system 
size at various temperatures for the 2D LL-model. (Monte 
Carlo simulation with 10® equilibrium steps (measured as the 
number of flipped Wolff clusters) after removing the first 10® 
MCS to achieve proper thermalization). 
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FIG. 7: Behaviour of the magnetic susceptibility versus the 
system size at various temperatures for the 2D LL-model. 
(Monte Carlo simulation with 10® equilibrium steps (mea- 
sured as the number of flipped Wolff clusters) after removing 
the first 10® MCS to achieve proper thermalization). 



D. Conclusion from the study of the m and x 
behaviours 

From this preliminary analysis, one could conclude 
that the 2D LL-model behaves similarly to the 2D XY- 
model. Namely, one can define a BKT temperature, T*, 
at which the system is critical, and a QLRO phase seem 
to occur at temperatures smaller than T*. 

We shall see below that this anticipated conclusion is 
refuted by more refined analysis. 



IV. CALCULATIONS OF THE CORRELATION 
FUNCTION EXPONENT 

In principle, very precise information can be obtained 
from the estimates of the fundamental critical exponent 
T] for the pair correlation function. This exponent may 
depend on the effective temperature T in the case of a 
continuous line of critical points. 

Several methods can be used independently in the case 
of the BKT transition. We shall use the following ones 
below: 

• Finite size scaling behaviour of the magnetization. 

At the critical point of a second order phase transi- 
tion, the magnetization m behaves as a power law 
of the system size, namely: m ^ L~^/'^ . This com- 
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bination of critical exponents is related to the ex- 
ponent which describes the critical decay of appro- 
priate pair correlation function, G{r) r^id-^+ri) ^ 
2(3 1 V — d—2+r], thus in two dimensions m ~ L^^^, 
Here, 77 is a function of T and FSS is used to es- 
timate the critical exponents at different temper- 
atures (in the low temperature phase) from figure 
© or figure dS]) ; 



• FSS behaviour of the magnetic susceptibility. 

The magnetic susceptibility, Xj scales as a power 
law involving the critical exponents 7 
namely:^ v 



relation [51 



the critical 
^ LT/''. Then, 



the 



and V, 
hyperscaling 



(3) 



holds, which gives an alternative estimate of the 
exponent 77. Here too, FSS may be used to extract 
an appropriate exponent from the data of figure ([3]) 
or figure ([7]); 

• Conformal Transformation (CT) method. 

Another tool can be used to estimate the value of 77 
for very large systems: the conformal transforma- 
tion method (see Cardy [1^). It is indeed a pow- 
erful method which applies to analyze any confor- 
mally covariant density profile (or correlation func- 
tion) from the mapping of infinite or semi-infinite 
geometries onto restricted geometries where simu- 
lations are actually performed [H, [s^l- We thus 
consider a density profile m(w) in a finite system 
with symmetry breaking fields along some edges in 
order to induce a non- vanishing local order param- 
eter in the 2D bulk of the system. Here, w stands 
for the location in the finite system. In the case 
of a square lattice of size L x L, with fixed bound- 
ary conditions along the four edges, using a simple 
mathematical transformation - called the Schwarz- 
Christoffel transformation [H, [s^ - one expects a 
power law in terms of an effective locally rescaled 
distance 



m{w) ~ ['^(w)] ^' 



(4) 



This method has been used extensively in [3a|40| to 
estimate the pair correlation exponent as a function 
of temperature for the 2D XY-model; 

• Systems in an applied magnetic field. 

According to the scaling hypothesis, the magnetiza- 
tion of a general d-dimcnsional system at a critical 
temperature, must follow the scaling law: 



(5) 



where H is the applied magnetic field. The mag- 
netic exponent is given in terms of the magnetic 
field scaling dimension yh-, S = yh/{d — yh), and 



G is a universal function [4l|. In the thermody- 
namic limit L~^H~^/y'^ <^ 1, the system follows 
the singular behaviour m H^^^ . In the other 
limit L^^H~^/y'^ ^ 1, the magnetization follows 
m ~ i^"'^'*. As we know that m ^ L^**/^ for the 
BKT transition, one deduces in this case the hy- 
perscaling relation under the form: 

d~yh = ri/2. (6) 

which again can be used to estimate the value of rj. 

A. The 2D XY-model 

The FSS methods are used from figure and figure 
([S])- The results obtained by the CT method are from 
the reference [lol- We add here a few words about the 
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FIG. 8: XY-model in a magnetic field H. For each tempera- 
ture the system sizes are: L — 32, 64, 96, 128, and the mag- 
netic field strengths are: H = 0.1, 0.01, 0.05, 0.001, 0.005, 
0.0001, 0.0005. The bold straight lines represent the large 
X = HL^" limit for a WeibuU-like PDF (W) for each tem- 
perature. The number of Wolff steps for each T, L and H is 
2 X 10^ 

2D XY-model in a magnetic field. The exponent ry at 
the BKT transition is 1/4, then yh — 15/8, and S = 15 
has the same value as for the d — 2 Ising model. Using 
equation ([5]), an excellent collapse of Y = H/m^ versus 
X = HLy^ at Tbkt = 0.893 has been obtained in [ll]. 
The authors obtained in this work the additional result 
that Y tends to a constant value in the limit HL'*"^ 
00. The saturation of Y for large magnetic fields means 
a WeibuU-like fieldless probability distribution function 
(PDF) [43] for the magnetization. We extended in the 
present work this result to other temperatures. 

Indeed, for a BKT transition there is a line of critical 
points for any T < Tbkt, and the equation ([5]) must be 
satisfied in all this range of temperatures. Then, we per- 
formed numerical simulations of the 2D XY-model in an 
applied magnetic field 44] for temperatures below Tbkt- 
Figure ([5]) shows the excellent collapse of the data, as well 



7 



0.3 



0.2 



0.1 





spin wave solution 


o 


L = 100 C.T. 


□ 


L = 200 C.T. 


▼ 


-n/2 

m~ L 


A 


,2-11 

z-i 


♦ 


from 




0.2 



0.4 0.6 



0.8 



FIG. 9: Dependence of the correlation function critical ex- 
ponent with the temperature for the 2D XY-model. (Monte 
Carlo simulations with 10^ MCS/spin after removing the first 
10^ MCS to achieve proper thermalization) . The conformal 



transformation (CT) values are from the reference [40 
dashed line is the spin-wave result 77 = kBT/2nJ. 
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B. The 2D Heisenberg model 

We compare now the values measured for i] for a system 
(the 2D Heisenberg model) which does not exhibit any 
transition, at least for temperatures above T = 0.1 
(see figure HI)). Then, we can expect the behaviours of 77 
to be quite different in the two cases. 
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as the saturation of Y for several magnetic field strengths 
and several system sizes, at all the temperatures < Tbkt 
investigated. The collapse of the data gives an estimate 
of the critical exponent 77 (T) using the hyperscaling rela- 
tion dH). 

This result invalidates a claim by Bramwell et al. [i^ , 
about the possible Gumbel-like distributions of the 
magnetization in the 2D XY-model at low temperatures. 
Indeed, the Gumbcl PDF does not lead to the saturation 
of the function Y in the limit HL^'^ 00. 

The figure ([9]) shows the values of the critical exponent 
77 for the 2D XY-model, calculated from FSS using the 
magnetization or the magnetic susceptibility, with the 
CT method (from [13]), and using the collapse of the data 
for the system in a magnetic field. All the methods lead 
to determinations of 77 which are in excellent agreement 
with each other. 

However, one can note that the behaviour obtained 
with CT is the most precise, because this method is able 
to reach the behaviour of the infinite-size systems, with 
system sizes computationally accessible, all shape effects 
being encoded (in the continuum limit) in the conformal 
mapping (hence the name of Finite Shape Scaling that 
we tried to introduce in Ref. [23]). The differences be- 
tween the results obtained from m and from x, are of the 
order of 5%. In the references [s^l and [i^, the authors 
used more than 6 x 10^ independent realizations for the 
estimate of the exponent rj. Realization of the hyperscal- 
ing relation ^ has been checked with a relative error 
smaller than 0.4%. 

More than simple illustrations, these results allow us 
to compare the cffeciency of the different methods for 
the proper definition of a quasi-long-range-order phase 
in systems with continuous symmetry. 



FIG. 10: Dependence of the effective correlation function crit- 
ical exponent with the temperature for the 2D Heisenberg 
model. The dashed line is the spin-wave result rj = ksT /nj 

The figure (fT0| shows the values of 77 calculated us- 
ing the FSS and the CT methods for the 2D Heisenberg 
model. The differences with the 2D XY-model are clear. 
It is particularly noticeable that the estimates given by 
the three methods are completely different. Even the 
shapes of the curves are different, and no agreement can 
be found. In particular, the determinations of 77 deduced 
from 777, and x, respectively, do not agree with each other. 
It follows that the hyperscaling relation ([3]) is definitively 
not satisfied. 

These observations give a strong evidence for the lack 
of any critical behaviour in the 2D Heisenberg model 
within the range of temperatures considered here. 



C. The 2D LL-model 

We present now a similar systematic work, performed 
on the 2D LL-model. In particular, the figure pT|) shows 
results for the 2D LL-model in an applied magnetic field 
at the expected critical temperature, T* . 

The results are consistent with the equation ([5]), and 
with the saturation of the quantity Y . The best collapse 
was obtained for 5 = 10.83. Then, using the hyperscaling 
relation ^ we found: 77 — 0.338 at this temperature. 
This result is shown in figure (jl2p and the agreement 
with the value of 77 obtained from the other methods is 
excellent at this temperature. 

In figure IT^ , we show the values of 77 for the 2D LL- 
model at various temperatures. Estimates were done us- 
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10000 



FIG. 11: 2D Lebwohl-Lasher model in a magnetic field {H = 
0.1, 0.01, 0.05, 0.001, 0.005, 0.0001, 0.0005) at T = 0.513. 
The bold straight line represents the large X — HV'^ limit 
for a WeibuU-like PDF (W). Each point corresponds to an 
average over 100, 000 independent realizations. 
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FIG. 12: Dependence of the critical exponent rj with the tem- 
perature for the 2D Lebwohl-Lasher model. (Monte Carlo 
simulations with 10^ MCS/spin after removing the first 10^ 
MCS to achieve proper thermalization) . The conformal trans- 
formation (CT) data are from reference [2^. The dashed line 
is the spin-wave result r] = AkaT /ZttJ . 



ing FSS of the magnetization (from reference [20| and 
figure O and of the magnetic susceptibility (from figure 
[7]). They are compared to the values estimated from the 
CT method (from reference I25j), and from the collapse 
of the data in a magnetic field. 

Nevertheless, one should stress that the agreement be- 
tween the values estimated by all the four methods is 
good, but imperfect. Small discrepancies can easily be 
noticed. In particular, the estimates which are based on 
the FSS behaviour of the magnetic susceptibility are sys- 
tematically smaller than those following from the other 
methods. The errors are larger than 10%, a relevant 
quantity probably far out of the range of acceptable er- 
ror bars. In a previous work [33 |. this error was only 
of 3%, but the authors then checked systematically the 



hyperscaling relation, equation ([3]), and they concluded 
that the hyperscaling relation is definitively not satisfied 
at low temperature (e.g. at T = 0.4) (see discussion in 
the section |V] below) . This has been the first evidence 
for absence of a QLRO phase in the 2D LL-model. 



V. REVISITING THE BEHAVIOUR OF THE 
2D LL-MODEL AT THE 'CRITICAL' 
TEMPERATURE 

In the previous section, we observed that the critical 
behaviour of the 2D XY-model and the behaviour of the 
2D LL-model were quite similar. Nevertheless, we saw 
small discrepancies for the LL-model, in particular in the 
behaviour of the second moment of the magnetization, 
and on the validity of the hyperscaling relation ([3|) . 

In reference [s^ . the authors performed intensive sim- 
ulations of the 2D LL-model at the temperature T = 0.4 
well below the expected critical temperature. Clear fail- 
ure of the hyperscaling relation, failure of collapse of the 
PDF in the first scaling form [s^, and analysis of the 
Binder cumulant, led to the conclusion that no QLRO 
phase exists in the 2D LL-model at this temperature. 

Here, we consider the problem of the existence of a 
critical at the temperature T* = 0.513, as this tempera- 
ture was reported as a possible BKT temperature for the 
2D LL-model (see section lllll and references (l^,!!!!). To 
elaborate on this question, we will report results about 
the probability distribution function, the magnetization 
scaling and the hyperscaling relation ([3]) at this tempera- 
ture. We shall also compare the stiffness for the LL- and 
XY-models. 



A. The first-scaling law 

If the 2D Lebwohl-Lasher model is critical at the tem- 
perature T* = 0.513, self-similarity must be observed at 
this temperature because no finite characteristic length 
can exist at criticality. Asymptotic self-similarity results 
in the so-called first-scaling law for the order parameter, 
here the magnetization m, [33| : 



(to)P(to) = $t(zi), 



with zi = 



(m)' 



(7) 



with $T a scaling function which depends only on the 
actual temperature T. Equation ([7|) is sequel of the stan- 
dard finite-size scaling theory |46j , but it is highly advan- 
tageous that equation ([7]) does not require knowledge of 
the values of any critical exponent. Note that, under this 
form, the hyperscaling relation. 



(to) 
a 



constant term, 



(8) 



is automatically realized with a the standard deviation of 
the order parameter. Indeed, one has: oc x/L'^i a 
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FIG. 13: PDF of the magnetization for the 2D Lebwohl- 
Lasher model at the temperature T* = 0.513, plotted in the 
first-scahng form ([7]). The scaling law is not confirmed for L = 
32 up to 1/ = 256. Wolff's single-cluster algorithm was used 
[sj. Each data set corresponds to average over 6,000,000 
independent realizations. 



scales then with the system size as: a Lt/^^^-i for 2D 
systems. Therefore, the ratio (m) /a should be a constant 
whenever the hyperscaUng relation ^ is satisfied |32] . 

Figure (|13p shows the behaviour of the PDF when the 
system size increases, for the 2D LL-model at the tem- 
perature T* — 0.513. We can observe failure of the col- 
lapse between these curves, specially for the largest sys- 
tem sizes, where the distributions tend to deviate from 
each other when L increases. One must conclude that the 
temperature T* = 0.513 is not critical for this model. 

In two others papers published recently, one can find 
the same kind of curve for the 2D XY-model at T = 0.6 < 
Tbkt [11 and £it T = Tbkt The collapse of all the 
PDF is excellent in the two result consistent with 

the known criticality of the system at Tbkt and all the 
temperatures below this value. In the same reference 
[33 |. the 2D LL-model was studied similarly at T = 0.4. 
Failure of the collapse was clear in this case, leading to 
the conclusion that the system is not critical below the 
assumed transition temperature. 

In 2003, Mondal et al. [30*1 reported a possible transi- 
tion temperature for the 2D LL-model at the temperature 
T = 0.548, larger than the values reported by Kunz et al 
in 1992 |24|. We thus also performed numerical simula- 
tions of the 2D LL-model at T = 0.548. In figure ^ 
the data are plotted in the first-scaling form ([7]). The 
data are far from collapsing. It is actually similar to the 
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FIG. 14: PDF of the magnetization for the Lebwohl-Lasher 
model at T = 0.548, plotted in the first-scaling form ([7]). Each 
data set corresponds to average over 6, 000, 000 independent 
realizations. 



behaviour of a system far away from any critical point. 
15r 




FIG. 15: m/a vs. L"^ at T = Tbkt for the 2D XY-model 
(top) and at T = r* = 0.513 for the 2D Lebwohl-Lasher 
model (bottom). For the XY-model, the data were obtained 
from the reference [i^. For the LL-model, the blue line is a 
power law fit. 



B. The hyperscaling relation 

In this section, we check numerically the hyperscaling 
relation under the form ([8]). In figure (fT5|) . {ra)/a is 
plotted versus 1/L for the 2D XY-model and the 2D LL- 
model respectively, at the BKT temperatures reported 
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or expected for each system, namely: Tbkt = 0.893 for 
the XY-model and T* = 0.513 for the LL-model. 

For the XY-model, the ratio is seen to converge to 
a definite value (~ 14.8) at the thermodynamic limit 
\/L = 0, while for the LL-modcl, a power law is the 
best fit consistent with the data of {m)/(j versus i, and 
this fit leads to a vanishing value of the ratio at the ther- 
modynamic limit. We conclude that the hyperscaling re- 
lation is very likely violated in this case, for the obvious 
reason that the system is not critical at the temperature 
T* = 0.513. 



C. The stiffness modulus 

Appearance of stiffness is an important consequence of 
the spontaneous symmetry breakdown in a system with 
continuous symmetry. It is closely related to correlations 
in the transverse response function, any spatial variation 
of the order parameter, perpendicular to the direction of 
the ordering, increasing the energy of the system. For 
this reason, the system generates a restored strength in 
response to any attempt to create a new configuration: 
this is known as stiffness. 

In the case of the 2D XY-model, this is the helicity 
modulus and it is defined as follows. Let 



7^ = -/3J^S,-g„S, 



(9) 



be the modified Hamiltonian in which cja is a rotation of 
angle (f> of the system around any axis, say here the x- 
axis. The difference between the free energy F{(j)) calcu- 
lated with the modified Hamiltonian, and the free energy 
calculated with the original Hamiltonian, 



F{c^) - F(0) 



1 /Tre-'^'^s- 



varies quadratically with the rotation angle 
for (/) <C 1. The coefficient, usually written: 
helicity modulus: 

F{(j)) - F(0) = -T(/.2 + 0(04) 



and it takes the form 




E 



(10) 

), at least 
T, is the 

(11) 



(12) 



where Jx is the rotation operator around the x-axis. 

In figure (fTB|) the stiffness modulus for the 2D XY- 
model is shown versus the temperature for various sizes 
from L = 16 to L = 512. The behaviour is typical of an 
order parameter: below Tbkt — 0.893, the plots are al- 
most independent of the system size, and this is the sign 



of the rigid phase with a finite value of the stiffness con- 
stant, while for temperatures larger than Tbkt = 0.893, 




FIG. 16: Temperature behaviour of the stiffness modulus T 
for the 2D XY-model. (Monte Carlo simulations with 10® 
MCS/spin after removing the first 10^ MCS to achieve proper 
thermalization) . 
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FIG. 17: Temperature behaviour of the stiffness modulus 
T for the 2D Lebwohl-Lasher model. (Monte Carlo simu- 
lations with 10*^ MCS/spin after removing the first lO" MCS 
to achieve proper thermalization). 



the system's response becomes weaker, leading to the de- 
crease of the stiffness modulus. The high-temperature 
phase is a disordered phase with the vanishing stiffness. 
The change of behaviour in the stiffness modulus between 
the two phases exhibits the appearance of a topological 
phase transition, with quasi-long-range order, but finite 
rigidity, in the low temperature phase. 

For the 2D LL-model, the same procedure can be used 
with eventually the following expression for the stiffness: 
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3J^ 



\<ij> 




(13) 



<jj> 




1000 



FIG. 18: StiflFness T versus system size L for three diflFerent 
temperatures for the 2D XY-model. T(I/) goes to a finite 
value in the limit L — > oo for T < Tbkt- (Monte Carlo 
simulations with lO" MCS/spin after cancellation of 10^ for 
thermalization) . 
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FIG. 19: Behaviour of stiffness modulus versus the system 
size for the 2D Lebwohl-Lasher model. (Monte Carlo simula- 
tions with 10^ MCS/spin after removing the first 10^ MCS to 
achieve proper thermalization). 



Then, we have plotted T versus the temperature as 
it is shown in figure (|17p for various sizes from L — 25 
to L = 400. Again, a behaviour similar to an order- 
parameter is observed in this case. At any temperature 
below r* = 0.513, one clearly identifies finite values of 
the stiffness independently of the system size. Above T*, 
the values of the stiffness drop to values close to 0. 

However, beyond the similarities, one can remark that 
the stabilization to a finite value of the stiffness in the 
low-temperature phase, is not similar in both models. 
One can see this point, in particular analyzing the de- 
pendences with the system size. In figure ([TSl) . T is 
shown versus L for three different temperatures of the 
2D XY-model. Two of them below the BKT transition 
(T = 0.5 and 0.7), and the other temperature just at 
the BKT transition. For the lower temperatures T goes 
fast to an asymptotic constant value. At Tbkt, the stiff- 
ness modulus reaches its saturation value algebraically. 
Then, it is clear that for the XY-model a sharp jump of 
the stiffness modulus is observed right at the topologi- 
cal transition. The behaviour is completely different for 
the 2D LL-model. In figure the dependence of T 

with the system size L is shown for various temperatures 
below T*. The stiffness does not reach a finite value in 
the limit L oo. This behaviour is similar to the 2D 
Heisenberg model |47| and to the fully frustrated Heisen- 
berg model [11] . In figure pJ]) , the slopes of the curves 
depend on the value of the temperature. It appears to 
be similar to the results of Winttel et al [i^ . A detailed 
study of the stiffness for frustrated spin systems can be 



found in [4£ 



VI. DISCUSSION AND CONCLUSIONS 

Our initial goal was to give definite conclusions about 
the possible appearance of a critical topological transi- 
tion, for example of the Berczinskii-Kosterlitz-Thouless 
type, in the 2D Lebwhol-Lashcr model. 

Our results exhibit a scenario more complex than ex- 
pected. It can be summarized as follows: 

• The behaviour of the order parameter (the mag- 
netization) with the temperature and the system 
size, pointed out a quasi-long range order phase for 
temperatures below the value T* = 0.513. One 
has noticed an excellent agreement of the correla- 
tion exponent rj estimated using conformal trans- 
formation, scaling of the order parameter with L, 
including scaling with a magnetic field. 

• If we investigate a quantity related to the second 
moment of the order parameter, that is the sus- 
ceptibility, a qualitative agreement with a QLRO 
phase behaviour is observed. However, the value 
obtained for 77 from this quantity does not compare 
well with the estimates from the order parameter. 
Actually, the hyperscaling relation ^ is not satis- 
fied for T <T*. This is in clear contradiction with 
the possibility of a line of critical point below T*. 
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Then, one has to conclude that a BKT transition 
is not observed for this model. 

• Although not studied in the present paper, the be- 
haviour of the Binder cumulant (which depends on 
the fourth moment of the order parameter) , as re- 
ported in [s^l, fully supports the above conclusion: 
no universality is observed at any finite tempera- 
ture. Signs of QLRO just disappear if larger mo- 
ments of the order parameters are used. 

• As the conclusion seems to depend on the order 
of the moment of the magnetization, it looks nat- 
ural to use the complete probability distribution 
function of the order parameter to clarify the situ- 
ation. Unlike the 2D XY-model, the 2D LL-model 
does not show any first-scaling-law collapse as it 
would be expected for a critical system. The ab- 
sence of self-similarity at any temperature discards 
any sort of transition in this model, which thus ap- 
pears more similar to the Heisenberg model than 
to the XY-model eventually. 

• The stiffness modulus mixes a second and a fourth 
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